NICOLAAS GOVERT DE BRUIJN, THE ENCHANTER 
OF FRIABLE INTEGERS 



PIETER MOREE 
In memoriam: Nicolaas Govert ('Dick') de Bruijn (1918-2012) 

Abstract. N.G. de Bruijn carried out fundamental work on in- 
tegers having only small prime factors and the Dickman-de Bruijn 
function that arises on computing the density of those integers. In 
this he used his earlier work on linear functionals and differential- 
difference equations. We review the relevant work and also some 
later improvements by others. 



1. Introduction 

The number theoretical work of Nicolaas Govert ('Dick') de Bruijn 
(1918-2012) comes into two highly distinct flavours: combinatorial and 
analytical. In his combinatorial number theoretical work de Bruijn 
even did some cowork with one of the all-time greats in this area: Paul 
Erdos (6 joint papers!). Here we will only discuss de Bruijn's work in 
analytic number theory. This was done mainly in two periods: 1948- 
1953 and 1962-1966. In the second period de Bruijn revisited his earlier 
subjects. Some of this work is joint with Jacobus Hendricus ('Jack') 
van Lint (1932-2004) 0131]. 

In sieve theory (see, e.g., [HI [55] or for a very brief introduction 
Section [6] below) one is interested in estimating in terms of elementary 
functions the number of integers in a prescribed sequence having a 
prescribed factorization. Usually one is interested in the integers in 
the sequence that are primes and if one cannot handle those, integers 
that have only few distinct prime factors. De Bruijn's work in analytic 
number theory belongs mainly to sieve theory, but he is not restricting 
the number of prime factors of the integers that remain after sieving. 
This is a milder form of sieving and here usually quite sharp estimates 
can be obtained. As a basic example one can take the friable number 
counting function ty(x,y). Let P(n) denote the largest prime divisor 
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of an integer n ^ 2. Put -P(l) = 1. A number n is said to be y-friabl^ 
if P{n) ^ y. We let S(x, y) denote the set of integers 1 ^ n ^ x such 
that P(n) ^ y. The cardinality of S(x,y) is denoted by ^(a;,?/). 
In 1930, Dickmari gl] proved that 

ty(x x l / u ) 

(1) lim w^x j 

where the Dickman function (today often also called Dickman-de Bruijn 
function) p{u) is defined by 



(2) p{u) 



1 for < it < 1; 

« Jo p('« — t) dt for u > 1. 



Note that p{u) > 0, for if not, than because of the continuity of p{u) 
there is a smallest zero u > 1 and then substituting u in (T5]) we easily 
arrive at a contradiction. De Bruijn's contribution was to provide a 
very precise asymptotic estimate for p{u) and in later works he pro- 
vided much more precise estimates then hitherto known for y). As 
a consequence various authors call p(u) the Dickman-de Bruijn func- 
tion and some authors ^(x, y) the Dickman-de Bruijn function. Both 
in depth and originality de Bruijn's work goes way beyond everything 
done on the subject until then. Only starting in the 1980's, when 
Hildebrand and Tenenbaum started to work in this area, various re- 
sults of de Bruijn were greatly improved upon. Most introductions of 
y) related papers mention some de Bruijn papers. The techniques 
de Bruijn used in studying p{u) are now standard in studying similar 
functions that arise in sieve theory, see, e.g. [Ml Appendix B]. 



^Somc authors use y-smooth. Friable is an adjective meaning easily crumbled or 
broken. 

2 Karl Dickman (1861-1947) wrote his paper when he was 69 years old, after 
retiring from the Swedish insurance business world. 
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In the rest of this paper we first describe de Bruijn's work on p(u) 
and then that on ty(x,y). Before describing de Bruijn's work on p(u) 
more in detail, we will discuss the basic properties of p{u) and then 
very briefly the linear functional work of de Bruijn. 

There is much more to friable integers than de Bruijn's work and 
its direct follow-up. The beautiful survey of Granville [51] makes very 
clear how many connections there are with other problems in (mainly) 
number theory, but also in probability, cycle structure of permutations 
(see Section H~3"j) . etc. etc. Here we will not go into this rich tapestry. 

It is assumed that the reader is familiar with the Landau-Bachmann 
O-notation (see wikipedia or any introductory text on analytic num- 
ber theory). Instead of log log x we sometimes write log 2 x, instead of 
(logx) A , \og A x. 

This paper is the extended version of a paper [HI] (in Dutch) written 
on request for an issue of Nieuw Archief voor Wiskunde dedicated to 
the mathematics and memory of de Bruijn. For a discussion of the 
work done on friable integers before de Bruijn, see Norton |83j and 
Moree (821. 



(3) p(u) 



2. Dick de Bruijn and the Dickman-de Bruijn function 
2.1. Elementary properties of p(u). From (T5]) we see that 

p ( u ) = 

u 

for u > 1 and thus an alternative way of defining p(u) is 

1 farO<ii<l; 
1 _ju £ (t-V Ldt forM>L 

It follows that p(u) = 1 — logw for 1 ^ u ^ 2. For 2 ^ u ^ 3, p(u) 
can be expressed in terms of the dilogarithm. Put r = (1 + v5)/2. We 
have, see, e.g., Moree [70] . 

7T 2 

p ( T 2 ) = 1 - 21ogr + (logr) 2 « 0.1046. . . 

60 

The author does not know any other explicit pair (u, v) satisfying u > 2 
and p(u) = v. 

Various papers discuss how to obtain high accuracy numerical ap- 
proximations to p(u), see, e.g., [70j [71]. The interest of de Bruijn was 
in the asymptotical behaviour of p(u) and similar functions. So that 
will be our main focus from now on. 

Since p'{u) = —p{u — l)/u for u > 1 and p{u) > for u > 0, it 
follows that p'{u) < and hence p{u) is monotonically decreasing for 
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u > 1. By (J2j) we find that up{u) ^ p(u — 1), which on using induction 
leads to p(u) ^ 1/M' f° r u ^ 0. It follows that p(u) quickly tends 
to zero as u tends to infinity De Bruijn was the first to obtain the 
following much more precise estimate for p(u), for u ^ 3, 

{-«{^"^«-^^+ (O 2 )}}- 

Note that (jl]) implies that, as m -)• oo, 

1 /e + o(l)\« 

formulae that suffice for most purposes and are easier to remember. 

De Bruijn obtained (jl]) as a corollary to his highly non-elementary 
Theorem [21 to be discussed below. In Section 12.81 an elementary proof 
of (@J is given. 

2.2. The Laplace transform of p(u). In de Bruijn's work the func- 
tion (e x — l)/x and related functions play a prominent role. To un- 
derstand its origin it is a very useful exercise to compute the Laplace 
transform of p(u). This transform is defined by 

POO 

p( s ) = / p(u)e- us du. 
Jo 

Using that < p{u) ^ l/[w]! it is seen that this integral is absolutely 
convergent for all complex s, and thus defines an entire fucntion of s. 
Using (j2J) we then find that 

/oo pi poo pu . 

up{u)e~ us du = J ue~ us du + J [J p{t)dVje- us du 

= jH Pit) ( f +% e— du) dt = 
It follows that p(s) = Cexp{— Ein(s)}, where C is some constant and 



"1- e -t 
Em(.s) := / dt, s G C, 



is the complementary exponential integral. The value of C can be de- 
duced by comparing the behaviour of p(s) and exp{— Ein(s)} as s tends 
to infinity along the positive real axis. On the one hand we have, by 
partial integration, 



POO 

/ p'(u)e- us du = -1 + sp(s) 
Jo 
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and since the integral, in absolute value, is bounded by 1/s, it follows 
that lim^oo sp(s) = 1. By integration by parts we see that 

r 1 e -t _ i r°° e -t r°° 

/ dt + / — dt = e _t logtdt = r'(l) = -7, 

Jo t Jit Jo 

where 7 denotes the Euler-Mascheroni constant and T(z) the Gamma- 
function. Thus, with s tending to infinity along the positive real axis, 

f 1 1 - e~* f s e~* 

Ein(s) = / dt- — dt + log s 

Jo t Ji t 



dt- — dt + — dt + logs 



t Ji * J s t 

= 7 + logs + 0(e" s ). 

It follows that C = e 7 . Thus we proved (see also [77J pp. 210-211] or 
[SSI PP- 370-372]) the following result. 

Lemma 1. We have 

f s 1 — e~* 

(5) p(s)=exp{7-/ — - — dt}. 

Jo 1 

As an amusing consequence it follows that for any ^ 8 ^ 1, 
p(0) = e 7 



poo 

\ p(t)dt = 5 + J2( n + 6 )p( n + 6 )- 



On applying the inverse Laplace transform to both sides in (jSJ) we 
obtain, for any real 00, 

■y i-ao+ioo r s e z _ 1 > 

(6) p(u) = -— / exp < 7 + / dz + us \ ds, 

2m Ja -ioc 1 Jo Z ) 

De Bruijn seems to have been the first to consider the Laplace transform 
of p{u) (in [23]). Lemma [1] is due to him. 

2.3. De Bruijn's £ function. An important quantity in the asymp- 
totic study of p(u) is the function £(u). For any given u > 1, £(«) is 
defined as the unique positive solution of the transcendental equation 

e ? - 1 

(7) ~T = U - 

It is easy to show that £(«) is indeed unique. Put f(x) = (e x — l)/x. 
Note that the Taylor series of f(x) is 

2 s 

/(x) = 1 + 2! + 3! + 4! +••• 
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So f(x) tends to 1 as x tends to 0+ and is strictly increasing in x. This 
clearly shows that f(x) = u has a unique solution x > for each u > 1. 
From (|7J) we infer, that as w — > oo, we have 

(8) £ = log£ + log M + 0(^-). 

We want to find an expression for £ = for it large. Note that 
for u sufficiently large 1 < £ < 21ogw. Therefore log£ = O (log log u). 
On feeding this into the right hand side of (jHJ) we obtain that £ = 
log it + O (log log u), which on feeding again into the right hand side of 
(JED gives 

4 = log u + log log u + O ( — . 

lOgM / 

Iterating often enough we can get an error term 0(\og~ k u) for any 
fixed k > 0. 

By means of standard techniques of asymptotic analysis (see [301 
2.3]) it is possible to go beyond this process of iteration. Using these a 
convergent series for £(it) can be given (see [65]). 

We will also need some information on /"£(£) dt. By integration of 
a precise enough estimate for £(w) one obtains 

( 9) /*{(*)* = «{ log „ + io g2 „ - 1 + + o ( ( !^) ») } . 

J i y. logw V x YOgU / J 

As we will see the similarity between ([9]) and (jlj) is no coincidence (see 
Theorem [3]) . 

We leave it as an exercise to the reader to show that 

lim u^'{u) = 1. 

u— >co 

2.4. Ramanujan's hand... On June 1-5, 1987, a meeting was held 
at the University of Illinois to commemorate the centenary of Ra- 
manujan's birth and to survey the many areas of mathematics (and 
of physics) that have been profoundly influenced by his work. One of 
the speakers, R. W. Gosper, remarked in his lecture, "How can you love 
this man? He continually reaches his hand from his grave to snatch your 
theorems from you." A few striking examples are given in an article by 
Berndt |12j . One example involves Dickman's result ([1]). Ramanujan 
[86l p. 337] made the following claim (written at least ten years before 
Dickman's paper appeared!) that is equivalent to the following claim. 
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Claim 1. Define, for u ^ 0, 

* 1 ,...,t fe ^i 

Then, for u ^ 0, ([I]) holds with 

00 (—-\) k 

k=0 

For more details we refer to the book by Andrews and Berndt [5j 
Section 8.2]. The asymptotic behaviour of the Ik(u) has been studied 
by Soundararajan [95]. Related formulae for p(u) as iterated integral 
can be, e.g., found in Buchstab [33], Chowla and Vijayaraghavan [3D], 
Goncharov [51] (see also p. 22 of |8J) and Tenenbaum J96J (12) with 
A = 1/u and t = 1]. 

Ramanujan in his first letter (January 16th, 1913) to Hardy gave an 
asymptotic for \l/(x, 3). See Moree [H2] for further details. 

2.5. Mahler's partition problem. Mahler's partition problem is to 
find an asymptotic formula for the number p r (n) of representations of 
an integer n in the form n = hq + n±r + r^r 2 + . . ., r > 1 a prescribed 
integer, in nonnegative integers no, ni, ri2, ■ ■ ., when n is large. By 
considering the functional equation 

(10) F'(y) = F(y/r), 

K. Mahler [7U [72] obtained an approximate formula for the number 
p r {n). De Bruijn [18] determined the asymptotic behaviour much more 
precisely. A generalization of de Bruijn's results was established by 
Pennington [85J, with l,r, r 2 , . . . replaced by a more general sequence 
Ai, A 2 , • • • satisfying a certain growth condition. 

By setting y = r x and F(r x ) = G(x), Mahler's functional equation 
is transformed into a functional equation of the form 

(11) G'(x) =e ax+l3 G(x- 1), 

with (3 real. In [2"5] , de Bruijn studied, using saddle-point techniques, 
the asymptotic behavior of real solutions of the above equation as x 
tends to infinity in case a > and (3 is a complex number. 

2.6. The linear functional equations papers. De Bruijn's initial 
mathematical interest was in combinatorics. According to J. Kore- 
vaar, de Bruijn was really excited about Mahler's partition problem 
(de Bruijn attended lectures of Mahler on this problem). This problem 
leads one to consider linear functional equations and the author would 
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not be surprised if this problem is the ultimate source of de Bruijn's 
strong interest in these type of equations. Also he had a Ph.D. student 
who wrote his Ph.D. thesis on this area of research [10] . 

In [19], de Bruijn studies the equation x~ a f(x) +f(x) — f(x — 1) = 
and shows that depending on a the behaviour of f(x) is quite different. 
He discusses an application to Mahler's partition problem (see Section 
12. 5p . In a later paper [20], de Bruijn treats the more general equation 
w(x)f(x) + f(x) — f(x — 1) = 0. It is shown that if w(x) is not too 
small as x tends to infinity (a sufficient condition is w(x) ^ 1/logx 
as x — > oo), then every solution tends to a constant /(oo), and an 
estimate is obtained for the difference \f(x) — /(oo)|. 

In [22], de Bruijn is concerned with equations of the type 

f{x)= [ K(x,t)f(x-t)dt. 
Jo 

He is interested in imposing conditions on K which guarantee that 
every continuous solution f(x) will be convergent, i.e., linx^oo f(x) 
exists and is finite (in such the kernel K(x, t) is said to be 

'stabilizing'). As an application he establishes the following result. 

Theorem 1. If F is continuous and satisfies 

(12) xF(x) = [ F(x-t)dt (x > 1), 

Jo 

then there exists a constant C such that F(x) = {C + 0(x~ 1 ^ 2 )}p(x) . 



De Bruijn used Theorem [T] in order to prove Theorem [2j 

2.7. De Bruijn's asymptotic formula for p(u). In 1951, de Bruijn 
[23J proved the following theorem. 

Theorem 2. As u — > oo we have 

1 r f^ u) e s - 1 "I 
P( u ) >= exjM 7 - u£(u) + / ds\. 

V27TM L Jo S J 

Note that the integral on the right hand side equals — Ein(— £). We 
remark that 

(13) < + Ein(-£) = J\(t)dt, 

which can be seen on noting that 

A« + Ein(-O) = (« - + £(«) = u > 1, 
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by the definition of £, and 

lim « + Ein(-£)) = 0, 

since — > as u — > 1+. Theorem [2] together with (TT3"]) and (JSJ yield 
the estimate @ for p{u). 

Below an outline of the remarkable proof of de Bruijn of Theorem [2j 
1) Note that 



If r f* e s - 1 

Fi( u ) = ttztt / cxp <^ — uz + / ds\dz 

'w Jo s 



2m 



satisfies f lT2|) . where the contour W consists of {z : $lz ^ 0, = — n}; 
followed by {z : = 0, — 7r ^ 3,2 ^ +vr}; followed by {z : $lz ^ 
0, 32 = +tt}. 

2) Use the saddle point method to estimate Fi(u) for large u. 

3) Note that for any continuous G satisfying the adjoint equation 

(14) uG\u-l) = G{u)-G{u-l), u> 1, 

and any continuous F satisfying uF'{u — 1) = —F{u — 1), w > 1, the 
function 

rb 



(F, G) = [ F(u)G(u) du - bF(b)G(b 



is a constant independent of b > 0. 

4) Note that 

+ / [ exp <uz — ds > — dz, u > —1, 

■00 1/ S «/ ^ ^ 

satisfies (fl4|) . 

5) Estimate Gi(-u), m large, by the saddle point method. 

6) Use the estimates of steps 2 and 5 to show that (Fi, Gi) = 1. 

7) By considering appropriate limits, prove that (p, Gi) = me 1 . 

8) By Theorem [T] we then infer that p(u) ~ e 7 Fi(u). 

9) Combine the latter result with the estimate obtained in step 2 in 
order to complete the proof. □ 

Applying the saddle-point method directly to (JH]) it is possible to give 
a rather shorter proof of Theorem [21 see [9H pp. 374 3760. Indeed, 
in this way a rather sharper version of Theorem [2] can be obtained, see 
Smida [M]. It is not clear why de Bruijn did not choose to follow this 
easier approach. 



3 The first to remark this was Tenenbaum in 1990 in the French original of his 
book. 
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In 1982, Canfield [35] gave a reproof of Theorem |2] using combina- 
torial analysis. He defines two step functions Pn(u) and gjv(u)- The 
stepsize for each is iV -1 ; the value of pjq on [jiV -1 , (j+l)A^ _1 ) is A(N,j) 
for j ^ 0; the value of q N on ((j - 1)JV _1 , jN' 1 } is B(N,j) for j ^ 1, 
where 

jjA(N,j) = ±(A(N,j - 1) + . . . + A(N,j - iV)) 

(15) A(JV, 0) = A(N, 1) = . . . = A(JV, N-l) = l; 
and 

^S(iV, j) = ^(B(N,j - 1) + . . . + B(N,j — N + 1)) 

(16) B(JV, 1) = £(iV, 2) = . . . = B(iV, AT - 1) = 1. 

The idea is that the step functions pn and qN are discrete models of 
(EJ). The generating series for A(N,j) and B(N,j) are not difficult to 
find: 

oo oo 2 N 

J2 B (N + l,j)x j =xJ2 A ( N J) xj =xexp{x + y + ... + ^}. 

Canfield then continues by showing that, for all A" and u, 

B(N, [uN]) ^ p{u) ^ A(N, [uN]). 

Using the circle method he then shows there exists an integer-valued 
function u(u) such that 

B(v(u), [uv(u)\) ~ A(v(u), [uv(u)]) 7 1== exp(7 - < - Ein(-£)), 

thus giving another proof of Theorem [3 

Alladi [1] generalized de Bruijn's result as follows (see also [98j P- 
374]). 

Theorem 3. For u ^ 1 we have 

p{u) = (l + O(i)) ^/^exp ( 7 - ^" £(t) dt) • 

Since lim^oo u^'{u) = 1 and because of the identity f[T3~j) . Theorem 
131 has Theorem [2] as a corollary. Alladi writes that "this improvement 
is in fact obtained by iterating de Bruijn's method a second time, more 
carefully!". He starts by showing that in de Bruijn's Theorem [1] the 
error term can be substantially improved. 

Hensley [59] gave another proof of Theorem [3j An arithmetic proof 
of Theorem [3] via the function ty(x,y) is contained in Hildebrand and 



NICOLAAS GOVERT DE BRUIJN, THE ENCHANTER OF FRIABLE INTEGERS 



Tenenbaum |63j. Hildebrand [62] established an asymptotic estimate 
similar to the one given in Theorem [3] for every function in a one 
parameter family of differential-difference equations. Smida |94j and 
Xuan 1 104] gave (quite technical) further sharpenings of Theorem [3j 

De Bruijn in his proof of Theorem [2] apparently was the first to use 
an adjoint differential-difference equation. Its application to questions 
arising in sieve theory was developed in several papers by H. Iwaniec 
and others, in particular in Iwaniec 



2.8. An elementary proof of de Bruijn's estimate fl4]). Evertse 
et al. IH pp. 109-110] showed that for u ^ 1 

(17) exp ( - J £(*) dt) < p{u) < exp ( - J £(t) di 

The short and elementary proof of ffTT|) only uses that p'(u)/p(u) is 
non- decreasing for u > 1. There is also an elementary argument for 
that, see, e.g., [781 PP- 35]. Using the elementary estimate @ we infer 
from ([T7D that 

p(u)=exp(-^ £(t)dt + 0(log«)). 
On inserting in this, the proof of @ is completed. □ 

3. Preliminaries on ^(x,y) 

3.1. Counting primes. In order to study \l/(x, y) it is essential to 
have some understanding of the distribution of prime numbers. As 
usual we let n(x) = J2 p ^ x 1 denote the number of primes p ^ x. We 
let li(a;) = J* dt/logt denote the logarithmic integral. Hadamard and 
independently de la Vallee Poussin in 1896 established the celebrated 
prime number theorem, stating that 

(18) tt(x) ~ , x — > oo. 

log a; 

It follows from the proof of Hadamard and de la Vallee Poussin that 

(19) n(x)=1i(x)+0(E(x)), 
with 

E(x) = xe~ ClVT ^, 

and Ci some postive constant. The currently best error term, due to 
I.M. Vinogradov and Korobov in 1958, is 

(20) E(x) =xexp{-c 2 log 3/5 x(log 2 £)~ 1/5 }, 
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with C2 some positive constant (see [3E1 Ch. 3]). Important to keep in 
mind is that often 

(2i) E/w «fM di . 



p^x 



E.g., we expect that J2 p ^ x \ ~ J2 j^i ~ Joglog^- Indeed, it is a 
classical result that 

(22) V i = log log x + c 3 + O(-L), 

•"^ p logx 

with C3 a constant. 

Writing n(x) = \i(x) + i?(a;) we obtain, using the Stieltjes integral, 

£ /(p) = £ /(*)<**(*) = [ ]^ dt+ [ f(t) dE(t) > 

making ( |2~T1) more precise. 



3.2. The Buchstab functional equation for \I>(£,y). Historically 
the first functional equation to be derived for the \l/(x, y) function 
was the Buchstab functional equation [33] (a few authors call this the 
Buchstab-de Bruijn equation), 

(23) *{x,y)=*(x,z)- V tt(-p), 

y<p^z 

where 1 ^ y < z ^ x. In the special case when y = 1 and z > 1 we 
obtain 

tt(z,z) = 1 + V*(-,p). 

In order to derive the Buchstab functional equation we start by observ- 
ing that the difference \l/(x, z) — \l/(x, y) equals precisely the cardinality 
of the set of natural numbers ^ x having greatest prime factor in (y, z}. 
Let p be some prime number. Notice that the cardinality of the set of 
natural numbers ^ x having greatest prime factor p equals ^(x/p,p). 
Thus, the Buchstab functional equation results. It can be used to 
determine ^(x,y) successively in the regions y ^ x, \fx < y ^ x, 
x 1 / 3 ^ y < y/x. In terms of u := log xj logy this corresponds with 
u ^ 1, 1<m^2,2<m^3.... Note that u is very close to the 
quotient of the number of digits of x and the number of digits of y. 
For u ^ 1 we obviously have ^(x,y) = [x] and so in particular 
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*ff(x/p,p) = [x/p] for y/x ^ p ^ x. On taking z = x in the Buchstab 
functional equation it then follows, for 1 ^ u ^ 2, that 

(24) ¥(*,*) = [*]- & 

Proceeding in this way one gets an exact expression for ^>(x,y) in 
increasingly large u regions. Suppose we have computed ^(x,y) up 
to u ^ h. For a term ty(x/p,p) in the Buchstab functional equation 
the logarithm of the first argument divided by the logarithm of the 
second argument equals (log xj log p) — 1 and this is < u — 1. So once 
we can compute *$?(x, y) up to u ^ h we can compute *&(x, y) up to 
u ^ h + 1. With each step of this process the expression for *ff(x,y) 
becomes more complicated. The idea now is to invoke analytic number 
theory to 'smoothen' the sums. From ( l24"j) and ( fl8|) we obtain 

(25) *(*, y) = x - V - + 0(tt(x)) = x - V - + 0(-^). 

' » p logx 

This together with (|22|) yields, for all x ^ 2 and 1 < w ^ 2, 

^x 1 /") =x(l-logu)(l + 0( r ^— )). 

logx 

This suggests that ^(x, x l ^ u ) ~ xf(u) for some function / for every 
u > 0. Heuristically this can be readily justified. In the Buchstab 
functional equation we put y = x x l u and z = x 1 !'" with u > v > 1 
arbitrary real numbers. We then infer that 

«/(«)»*/(«)- £ f/(||f-i) 

x 1 / u <n^x 1 /" 

_ , v / x Aogx v dt 



l/« t log i log t 

r f( W - 1) ^ 

= xj(v)— x / aw, 

where in the second step we used the approximation principle pip and 
in the final step the transformation t = x x l w was made. Assuming that 
the involved error terms on dividing by x tend to zero as x and y tend 
to infinity, we conclude that the function f(u) satisfies 



f{u) = f(v) - / dw. 



w 



Thus f(u) = p(u). 
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3.3. The Hildebrand functional equation for ty{x,y). The Hilde- 
brand identity [60J reads 



V(x,y)logx = f ^^dt+ V V(-^,y)\ogp. 



p^y 



The advantage over the Buchstab identity is that one of the parameters 
is held fixed and that the coefficients are non-negative. It was the use of 
this identity that led to serious improvements of several of de Bruijn's 
results. For most purposes it is enough to work with the approximate 
identity 

V(x,y) logx = Y]*(-,y) logp + 0(E(x,y)), 
p^y 

where one can take E(x,y) = x. For all y ^ log 2+e x one can take 
E(x,y) = *ff(x,y). An heuristic proof of *ff(x,y) ~ xp{u) based on the 
Hildebrand functional equation can also be given, see [6U pp. 416-417]. 

3.4. Rankin's method. In his work in the 30's on the gaps between 
consecutive primes Rankin [87J introduced a simple idea to estimate 
*&(x, y) which turns out to be remarkably effective and can be used in 
similar situations. Starting point is the observation that for any a > 

(26) *(s,y)< E ^ = * r tt*,v), 

n£S(x,y) P(n)^y 

where 

C(s,y) = U(l- P -T\ 

p^y 

is the partial Euler product up to y for the Riemann zeta function ((s). 
Recall that, for 9Rs > 1, 



oo 

cw = E^ = IIr jL 



p 

n=l p 



Note that the logarithmic derivative of x s ((s, y) with respect to s equals 

logp 



log x - ^ 



p s - 1 

p^y 

It is minimized for a = a(x,y), where a(x, y) is the (unique) solution 
of 

> ; r = lOg X. 

p^y 
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From a sufficiently sharp form of the prime number theorem one, see 
[55] . can derive the estimate 

log(l +y/\ogx) ( /loglog(l + y) 

a{x, y) = 1 + \ 

logy V V logy 

The upper bound x a ((a, y) is only too large asymptotically by a factor 

a^2n(j) 2 (a,y), 

where <j>k(s, y) is defined as the fcth partial derivative of log ((s, y) with 
respect to s. This factor turns out to be relatively small. This follows 
from one of the deepest results in this area, due to Hildebrand and 
Tenenbaum [63J, namely that uniformly in the range x ^ y ^ 2 one 
has 

x a ((a 1 y) ( 1 ^/logjr 



»(«,»)= x :rJ ] ,(i + Q(-)+o(- 



Of course, in f[2T)j) one is free to make any choice of a. E.g., the choice 
a = l — l/(2 log y) leads to 

C(*,V) « exp { £ ^} < ex P { E E ^) } « lo ^' 

p^y p^y p^y 

which gives rise to ^(x, y) <C xe~"/ 2 logy. The factor logy can be 
easily removed, see [98j p. 359]. 

As a further example of working with Rankin's method, let us try to 
estimate $?(x, \og A x) for A > 1. Letting a = 1 — l/A, we get 

>/ \ rr sr^ P 1//A V l ^ A logx 

logC a,y < Vp- ff = E ^ T < 1 — f • 

^ ^ p logy loglogx 

We conclude that 

(27) ^Or,log A x) < ^i-iM+od/iogiogx)^ 

We will see shortly that the upper bound is actually sharp. 

3.5. A binomial lower bound for ^>(x,y). Let 2 = p\ < p 2 < P3 < 

. . . be the consecutive primes. Let pj~ be the largest prime p ^ y (thus 
k = 7r(y)). Evidently n is in S(x,y) if and only if we can write n in 
the form n = p e \p e 2 2 ■ ■ • Pk k > where the are non- negative integers and 
n < x, that is, 



(28) eilogpi + e 2 logp 2 H h e fc logp fc ^ log a;. 
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Since each log pi ^ logy, any solution to e\ + e2~\ — - + ek ^ [j^ff] = N 
gives a solution to ( 128]) with x = y u . It follows that 



(29) 



e 1 + ... + e k Hiu] \ ^' / 



(If the reader cannot prove the last equality, (s)he can consult [77J p. 
208].) 

4. De Bruijn's work on ^(x,y) 

4.1. De Bruijn's A function. Let us define (as many authors do) 
p{u) = for u < 0. In [23] de Bruijn introduced the function A(x,y). 
He defines it as follows: 

(30) Hx , s) = x r p (^^n) d{ ^ 



jo ' log?/ / t 

Since p(tt) = for u < and p(w) is continuous everywhere apart 
from the jump at u — 0, the integral is well-defined, unless the jump 
of p coincides with a jump of [t]/t. This is the case if x is a positive 
integer n; this turns out to be a minor nuisance and can be dealt with 
by defining A(n,y) = A(n + 0,y). The function A(x,y) satisfies an 
equation which is similar to the Buchstab equation (123"]) . It arises from 
( 12"3"|) by writing the sum in that equation formally as a Stieltjes integral 
J ty(x/t,t)d,7r(t) and replacing n(t) by li(£). One can check that indeed 

f z x dt 

(31) A(x,y)=A(x,z)- / A(-,t)- Q 



t logt 



Further, for < x ^ y we have 

A(x,y) = x d(^) =x v -^ = [x] =V(x,y). 
Jo t x 

Thus A(x, y) satisfies the same initial condition as ^(x, y) and obeys a 
functional equation which is the friable version of Buchstab's identity. 

De Bruijn showed that A(x, y) closely approximates *&(x, y) in a cer- 
tain (x, j/)-region (this region was substantially extended in 1989 by 
Saias [SD])- Further he shows that both A(x,y) and $f(x, y) satisfy an 
expansion of the form 

X / , a rl r h E n [X), 

' log y 

r=0 to y 

with E n (x) a term he describes rather explicitly and 

K(z) = CLq + d\Z + aiz 2 + . . . 
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the Taylor series around z = of z((l + z)/(l + z). The term E n (x) 
is only an error term if u is not close to an integer from above. For a 
deeper understanding of this phenomenon, see Theorem 1.1 and Lemma 
3.1 in Hanrot et al. [57] . 

4.2. De Bruijn's uniform version of Dickman's result. In [24] de 
Bruijn gives estimates for the differences 

y) — A(zc, y)\ and \A(x, y) — xp(u)\. 

These estimates involve the error term E(x) in the prime number theo- 
rem, cf. ffl~9]) . These bounds, when combined with the sharpest known 
form of the prime number theorem (having error term (|20l)). yield the 
following uniform estimate of Dickman's limit result ([Q). 

Theorem 4. The estimate 

(32) * {x , y)=xM {l + Op&±Q)}, 

holds forl^u^ log 3/5 ~ € y, that is, y > exp(log 5/8+e x). 

Hildebrand [61] improved this substantially to the range 

1 ^ u ^ exp(log 3/5_e ?/), that is, y > exp((loglogx) 5/3+e ). 

There are limitations as to how far this range can be made. Hildebrand 
[BTi] showed that the estimate (13"2"|) holds uniformly for 

1 «C u sC y 1/2 ~% that is, y ^ log 2+e x, 

if and only if the Riemann Hypothesis is true. 

One might wonder in which region $?(x,y) ^> xp(u) is valid, as did 
Hensley [58J . Canfield et al. |36j proved that there is a constant C such 
that ^(x, y) ^ xp(u) exp{C(log 2 uj logw) 2 } uniformly for x > 1 and 
u > 3. It seems that "$?(x, y) ^ xp(u) for a very large (x, y)-region. 

4.3. The average largest prime factor of n. Recall that P(n) de- 
notes the largest prime factor of n. De Bruijn [24] applies his results 
on A(x, y) to show that 

(33) VlogP(n) = Xx\ogx + 0(x), with A = / . p(yU \ du. 
V Jo (1+m) 2 

By partial integration it follows that 

limi y lQ f P(n) =A, 
x-xx x ^— ' log n 

thus proving an heuristic claim by Dickman [41]. The interpretation 
of the result is that for an average integer with m digits, its greatest 
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prime factor has about Am digits. De Bruijn's estimate can be further 
sharpened to 



see [981 Exercise III. 5. 3]. The constant A is now known as the Golomb- 
Dickman constant. Mitchell [76j computed that A = 0.62432998854 .... 
Golomb et al. [5U] have defined a constant p which is related to the 
limiting properties of random permutations. Let L n be the expected 
length of the longest cycle of a random permutation of n letters. Define 
p n = L n /n. (Thus n\ — 1, p.% — 3/4, /i 3 = 13/18, p^ = 67/96.) It 
can be shown that the numbers p n are monotonically decreasing with 
n and thus a limit p exists. It turns out that A = p. Thus on average 
the longest cycle of a permutation of n elements has length An. Much 
earlier Goncharov announced in [51] and proved in [52] that the ratio 
of permutations in S n having longest cycle ^ n/u, with u > 1 and 
fixed, tends to p(u) as n tends to infinity. Goncharov, unaware of the 
occurrence of p{u) in number theory, also expressed p{u) as an iterated 
integral, see Arratia et al. [8, p. 22]. Shepp and Lloyd [93] derived a 
second formula for A: 



It is now known that the cycle lengths in a permutation are typically 
Poisson [53] . 

It is striking that p(u) turns up both in the setting of greatest cycle 
length and greatest prime factor. It also shows up in the setting of 
greatest irreducible polynomial factor of a polynomial over a finite field, 
see, e.g., Car [37]. For a unified treatment see Arratia et al. [7J (or 
their book [S]). 

Meanwhile there are many results (see, e.g., [131 E2]) involving P(n), 
one of the more interesting, see [13], being 



Later Scourfield [92] improved on the error term. Sums of the form 

En<i/(-P( n ))i where / is some arithmetic function can be evaluated 
using the identity 



together with sufficiently sharp estimates for ^(x,y). 

A quite general result is due to Tenenbaum and Wu [101] . who gave 



logP(n) = Ax logo; — A(l — 7)0; + 0{x exp(— log 3 / 8 e x)), 
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an asymptotic expansion for ^2 n<x f (n) {log P(n)} r , where / is a fairly 
general non-negative arithmetical function (which need not be multi- 
plicative) and r is any positive real number. 

4.4. De Bruijn's 1966 estimate for log\I/(x,y). We now turn our 
attention to de Bruijn's 1966 paper [27]. In it he is concerned with 
estimating log ^f(x, y). Set 

(ia\ v logs / y \ y ( logx\ 

34 Z := log 1 + - +- log 1 + . 

logy V logx/ logy V y J 

Theorem 5. We have, uniformly for x ^ y ^ 2, 

\og^(x,y) = Z{ 1 +o(^- + -+ 1 - )}■ 
L Vlog y u log log 2x / J 

Actually, in the error term u~ 1 can be left out, see Tenenbaum [98, 
pp. 359-362] for a proof. Let us apply Theorem [5] with y = log^x and 
A > 1 fixed. One then finds that 

V(x,\og A x) =x i-iM+o(i/iogio g ,) ; 

thus improving on (j27J. Already the results of de Bruijn in [23] imply 
as a special case that \l/(x, \og A x) = x 1 ~ 1 / j4+ °( e ). Following de Bruijn's 
argument very closely Chowla and Briggs [5H] gave a reproof that they 
believed to be considerably simpler. 

The proof of Theorem is very elementary (but takes several pages 
to work out). The upper bound is deduced from Rankin's upper bound 
( 126]) with an optimal choice of the parameter a . The lower bound is 
based from the inequality ( 129]) together with a lower bound for the 
right hand side of ( 129]) obtained by Stirling's formula. 

Theorem [5] clearly shows that there is a change of behaviour of 
^(x,y) as y ~ logx. If y/\ogx —> oo, then the first term in (134"]) 
dominates, whereas Z is asymptotic to the second term in ( I34p when 
y = o(logx). The change in behaviour is due to the fact that if y is 
small compared to logx, then many prime factors of a 'typical' integer 
in S(x, y) occur to high powers, whereas for larger values of y most 
prime factors occur only to the first power. P. Erdos [32] showed that, 
if x —>■ oo, we have 

i ,w i \ (log 4) logx 

log M/(x, logx) ~ — — • 

log log x 

5. De Bruijn and the sieve of Eratosthenes 

Define $(x, y) as the number of integers n ^ x having no prime 
factors ^ y. In the sieve of Eratosthenes all the integers 1 ^ n ^ x, 
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which are multiples of the primes 2 ^ p < y/x are removed. What 
remains after this process is the number 1 and all the prime numbers 
p in the range \fx ^ p ^ x. It follows that 

(35) $(x, y/x) = 1 + 7r(x) - ir(y/x). 

Buchstab 1 33 1 showed that 



(36) &(x,y) ~xe~<uj{u) JJ(1 - - 



p 

p<y 



where u(u) = 1/u for 1 ^ u ^ 2 and, for all u ^ 2, 



uuj{u) = 1+1 u(t)dt. 



u-l 



By a theorem of Mertens we have n p <y(l — Vp) 1 ~ me 7 / l°gy> an d 
thus we can reformulate (1361) as 

x 

$(x, y) ~ Ubj(ll) 



logx 

It follows, e.g., that $(x, a/x) ~ x/logx, an asymptotic one also finds 
from ([HD and (155} . 
De Bruijn [21] writes 

$(x, y) = x JJ(1 - jT 1 )^^, y) 

and shows that 

y*- M w(t)dt| < ci exp(-c 2V / logy), u > 1, j/ ^ 2, 

for appropriate constants c\ and c 2 . 

De Bruijn [201 Example 1] showed that lim^oo u{u) exists, and, if 
we denote it by A, then 

(37) w(u)=A + 0(-L). 



In a later paper de Bruijn [21] showed that A = e 7 . De Bruijn stated 
that the right order of magnitude of \u(u)— e _7 | "is probably something 
like 

(38) exp( — u log u — u log 2 it) " . 

This was proved subsequently by L.K. Hua [66] by an "advanced cal- 
culus argument" and by Buchstab [34| using an arithmetic argument. 
The formula (138]) reminds one of the asymptotic estimate (j3j) and in- 
deed it turns out to be natural to write uj{u) — e~ 7 = p{u)h{u) in order 
to give a sharp estimate for the difference, see Tenenbaum [99j Lemme 
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4]. A further connection between u(u) and p(u) is that their Laplace 
transforms p(s) and u(s) are related by 

sp(s) 

A surprising and important application of the properties of a; (it) has 
been mady by Maier [73]. By an ingenious construction, Maier proved 
that the number of primes in short intervals [x,x + log a?], C > 1, 
is sometimes larger than, and sometimes smaller than the expected 
number log* 7-1 x. Put W(u) = e 7 c<;(it) — 1. Put M + (v) = max u ^ v W{u) 
and M-(v) = mm u ^ v W(u). Maier proved for any fixed C > 1, 

r ir(x + \og c x) -n(x) 

hm sup ^— : ^ 1 + M_|_(G); 

2-s.oo log x 

„ tc(x + log c x) — tc(x) , , 

lim inf K —^- ^ 1 + M_(C); 

x^oo log x 

Using a method due to de Bruijn involving the adjoint equation of 
a; (it), Maier showed that W(u) changes sign in every interval of length 
one. Hence, for all C ^ 1, M + (C) > and M_(C) < 0. The paper 
of Maier led to a lot of follow-up work on the lack of equi-distribution 
of the primes. In this work methods from the study of friable integers 
and differential- difference equations play an important role, see e.g. 

6. HOW SPECIAL ARE p AND CU? 

The reader might wonder how special p(u) and u(u) are. The an- 
swer is that many similar functions occur in number theory, especially 
in sieve theory (the sieve of Eratosthenes being the easiest example of 
a sieve). Let us give some examples. 

1) Let Q be a set of primes, and denote by ^(x,y,Q) the number of 
positive integers not exceeding x that have no prime factors from Q 
exceeding y. Suppose that the number of primes 7r(x, Q) of primes 
p ^ x in Q satisfies 

ti(x,Q) = 6\i(x) + 0(—^), 

log X 

with < 5 < 1. Then Goldston and McCurley showed that 
*(x 7 y,Q)=xT S (u)(l + 0(^-)) 1 
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uniformly for u ^ 1 and y ^ 1-5. Here 75 (tt) is a function similar to 
it is the unique solution of Ts(u) — 1 for ^ u ^ 1 and 

ut' 5 {u) = —5ts(u — 1) for w > 1. 

This family of differential-difference equations had been earlier inves- 
tigated Beenakker [10J, a Ph.D. student of de Bruijn. 

The latter result is not valid if Q is the set of all primes (then 5 = 1), 
as the log(w + 1) in the error term in (132]) is sharp and cannot be re- 
moved. 



2) Let k ^ 1 be an integer. Let Pk{n) denote the kth largest prime 
factor of n. Let ^(x, y) denote the number of integers n ^ x such 
that Pk{n) ^ y. Knuth and Trabb Pardo [69] showed that 



hm = pi,k{u), 

x—toc X 

where p\ t k{ u ) is a function similar to the Dickman-de Bruijn function. 

3) In sieving we consider a finite sequence A of integers and a set P 
of primes. Let X denote the number of integers in A. The goal is to 
count the number of elements S(A, P, z) that remain after sifting A by 
all the primes in P less than z, that is we are interested in estimating 

S(A,P,z) :=\{aeA:(a,P(z)) = l}\, 

where P(z) := ripep P < Z P- This formulation is much too general if 
one wants to obtain non-trivial results and so we impose regularity 
conditions for d\P(z) on the subsequence 

A d := {a G A : a = 0(mod d)}. 

We assume that we have approximations of the form 

(39) \A d \ = ^X + r d , 

d 

where cu(p) is a multiplicative function satisfying < cu(p) < p such 
that r d is a remainder that is small. For p (jL P we put oj(p) = 0. We 
define 

^w-n(i-^)- 

Intuitively V(z) can be regarded as the probability that an element of 
A is not divisible by any prime p < z with p in P. We further assume 
there is a positive constant c\ such that 

(40) < ^ < 1 - 1 

p Ci 
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for every p in P and that there are constants k > and c 2 > 1 such 
that 



< c 2 , for all w 2 > W\ > 2. 



41 I V -^bgp-Klog(— 
I *— ' p W\ 

W\<p<W2 

The constant k is called the dimension of the sieve and is roughly the 
average of u (p) over the primes as 

logp 



~ log w 2 



(43) S(A, P, z) > XV(z){f K (-M) + 0} + R y , 



p 

p<W2 

as W2 tends to infinity. Under the assumptions ( |39l) . ( )40l) and f j411 one 
obtains the following inequalities for S(A,P,z): 

(42) S(A, P, z) < XV(z){F K (^-) + 0} + R v , 

log z 

and 

dog z' 

respectively, where F K {u) and f K {u) are sieve auxiliary functions de- 
pending on the dimension k, y is a free parameter, o is an error term 
and R y is an error term formed of the r^'s of (I3"5|) . cf. the inequality 
(|44p . In applications the free parameter y is chosen in such a way that 
the remainder term R y is small in comparison with the main term con- 
taining X. 

The Buchstab functional equation (123]) has a counterpart in sieve 
theory: 

S(A,P,z) = S(A,P, Zl ) - Yl S(A,,P,p), z>zx>2. 

zi<p<z, pSP 

From this we can obtain a lower bound for S(A, P, z) on assuming that 
a lower bound for S(A,P, zi) is already known, and that we have an 
upper bound for the terms in the sum. For z\ small enough an estimate 
for S(A, P, z) follows from the so called Fundamental Lemma of Sieve 
theory which states that under the regularity conditions (|39|) . (140]) and 
(|4T]) supplemented with the condition that \r d \ < oj(d) for d\P(z) and 
d square free, we have for z\ < X, 

S(A, P, z x ) = XV{ Zl ){l + 0(e-"/ 2 )}, 

where u = logX/log^i [551 Theorem 2.5]. 

As a demonstration of the above let us consider the Selberg's upper 
bound sieve (55J [89]. Let k > 0, u = logy/log^ > and u(d) the 
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number of distinct prime divisors of d. Under the conditions ( 1391) . 
(HDD, (JHJ), we have 

a K (u) V logy 7 1 d< ^ m 

where <J K (u) is the Ankeny-Onishi-Selberg function. The implied con- 
stant depends only on k, c\ and c 2 . One has a K (u) = p K {u)u K ~ x , where 
Pk(u) = 1 for < u < 1 and 

p' K (w) = - u)"~V(w -!),«> 1- 

Note that o"i(tt) = p(w). It seems that cr K (w) was introduced indepen- 
dently by Ankeny and Onishi [5] and de Bruijn and van Lint [25] . 
The function cr w (u) satisfies the differential-difference equation 

(45) uf(u) + af(u) + bf(u-l) = 0, 

with a = 1 — k and b = k. The Buchstab sifting function w(-u) satisfies 
(T4"5"|) with a = 1 and 6 = — 1. Solutions to the equation ([415]) appear 
to have been first studied by Iwaniec [57] in connection with his work 
on Rosser's sieve. A more systematic investigation was carried out by 
Wheeler |1U3] . Both Iwaniec and Wheeler make intensive use of the so 
called adjoint equation 

(46) ug\u) + (1 - a)g{u) - bg{u + 1) = 0, 

which, in some sense, is easier to deal with than (f4"5]) . Whereas the 
functions satisfying (1451) can exhibit a rather erratic behaviour, it turns 
out that there exists a solution of (1461) which is analytic in the right half 
plane. Since solutions to the equations (j4"5]) and (|46l) are connected to 
a simple integral relation (the right hand side of (T4T1) . which turns out 
to be a constant), one can derive asymptotic information for a given 
solution to (|45l) by studying the asymptotic behaviour of a suitable 
solution to (|46l) . Notice that this is precisely the approach followed by 
de Bruijn in proving Theorem [2] It has been established by Hildebrand 
and Tenenbaum [05] that also in this general case a much more direct 
approach is possible (cf. the paragraph in Section 12.71 following the 
proof of de Bruijn). 

The problem of describing the general solution to (|45l) amounts to 
describing solutions f{u) = f(u;(p) satisfying (T4"5]) for u > uq and 
f{u) = (p{u) for m — 1 ^ m ^ m , where uq is any positive real number 
and 4>{u) is any given continuous function on [w — l,«o]- Given two 
functions / and g defined on [u — l,u ] and [u ,Uq + 1], respectively, 
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we set 

(47) (/, g) = u f(u )g(u ) - b / f(u)g(u + 1) du. 

JUQ — l 

If / and g are solutions of ( 145]) . respectively (|46|) . then this "scalar 
product" is independent of u . 

Hildebrand and Tenenbaum [65J showed that there is a solution 
F(u;a,b) and that there are solutions F n (u;a,b) for every integer n 
that form a basis of the solution space, that is any solution can be 
expressed in the form 

oo 

f(u) = aF(u;a,b) + ^ a n F n (u;a,b), 

n=— oo 

with suitable coefficients a and a n depending on the initial function 
cf). De Bruijn [25] by a different approach had earlier found a similar 
result for solutions of the differential-difference equation (|TTT) . 

Earlier special cases of f|45|) had been considered by Alladi [I] , Beenakker 
[TO] , Hensley [59], Hildebrand [62] and Wheeler |103j amongst others. 

7. Arithmetic sums over friable integers 
In [29] de Bruijn and van Lint consider sums of the form 

ng5(x,y) 

with / a non-negative multiplicative function. Many of the usual tech- 
niques can be applied here, e.g., the analogues of the Buchstab and 
Hildebrand functional equation hold, see e.g. [80J. 

For an enlightening historical introduction and an extensive bibliog- 
raphy on the subject, see Tenenbaum and Wu [lOOj . For the state of 
the art of this area see Tenenbaum and Wu jlOlj . 

Let ^f m (x,y) denote the number of ?/-friables that are ^ x and co- 
prime to m. It turns out that evaluating the ratio ^ m (x/d,y)/^(x,y) 
forl^rf^x, m ^ 1, x ^ y ^ 2, is a crucial step for estimating arith- 
metic sums over friable integers, see de la Breteche and Tenenbaum 
[15]. Alladi [21 13] noticed a duality which shows that *$?(x, y) is related 
to the sum of the Mobius function over the uncancelled elements in 
the sieve of Eratosthenes, and the sum of the Mobius function over the 
smooth numbers is related to y). 

An interesting related theme is that of the friable Turdn-Kubilius 
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inequality. The Turan-Kubilius inequality 

J2\f(n)-A f (x)\ 2 <^xB f (x), 

holds uniformly for strongly additive functions /, where Af(x), Bf(x) 
stand for certain sums over primes p ^ x, depending on /. If 

B f (x) = o(A f (x), 

it can be shown that / has normal order Af(x), that is, for any e > 0, 
we have 

\f(n)-A f (n)\^e\A f (n)\ 
on a set of integers n of density 1. One can try to obtain an inequality 
analogous to the Turan-Kubilius inequality for the corresponding sum 
over the elements of the set S(x, y). The first result in this direction is 
due to Alladi PQ , who found that he needed a stronger version (Theorem 
12]) of de Bruijn's Theorem [5] for this purpose. For further developments 
the reader is referred to works of Tenenbaum and his collaborators 

dang mi Ea Eg. 

8. De Bruijn's analytic number theory work in a nutshell 

-HE], 1948- 

Mahler's partition problem is to find an asymptotic formula for p r (n), 
the number of partitions of n into powers of a fixed integer r > 1. De 
Bruijn gives an asymptotic for log p r (n) improving on Mahler's result 
HQ (cf. Section^. 

-[EH], 1949- 

This paper is mostly in pure analysis. At the end of the paper there is 
a brief discussion of connections with Mahler's partition problem and 
a promise to come back to this in a future paper (|25j). 

-[20], 1950- 

As an application of more general results de Bruijn shows that uj{u) = 
A + 0(l/[w]!) for some constant A. 

-j2Tj, 1950- 

He gives a sharp estimate for $>(x,y). Further he gives the first proof 
that the Buchstab sifting function u(u) tends to e -7 as u tends to 
infinity and formulates a conjecture (see ( 13 8 p that was later proved by 
'the father of modern Chinese mathematics' (Loo-Keng Hua) |66j . 
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-[22], 1950- 

He proves Theorem [TJ 

-[23], 1951 ~ 

The main result he proves here is Theorem [21 further Lemma [I] and 
the inverse Laplace transform formula (jUJ). As a corollary of his main 
theorem he obtains (j3J). Further, he introduces the function £(tt). 

-[21], 1951- 

He introduces A(x,y), shows that it closely approximates both *&(x, y) 
and xp(u). These estimates in combination with the error estimate (T2"U1) 
in the prime number theorem lead to a uniform version of Dickman's 
result ([1]), namely Theorem HI Further, he uses his A(x, y) function to 
prove the estimate fl33|) . 

[22], 1953- 

He studies, using saddle-point techniques, the asymptotic behavior in 
x of real solutions of the differential- difference equation 

F'(x) = e ax+ ^F{x - 1), with a > 0, G C. 

This is related to Mahler's partition problem (cf. Section 123]) . in which 
case (3 is real. 

-[26], 1962/1963- 

If n and x are positive integers, let f(n, x) denote the number of integers 
1 < m < x such that all prime factors of m divide n. Let j(k) = Yl P \kP 
denote the squarefree kernel of k. Erdos in a letter to de Bruijn conjec- 
tured that F(x) := J2 n ^x f( n ^ x ) satisfies log(F(a;)/x) = 0(log 1 ^ 2+e x). 
Put G(x) = J2 n <x f( n ' n )- ^ e Bruijn, using a Tauberian theorem by 
Hardy and Ramanujan, proves that 




thus proving a stronger form of Erdos conjecture. (Later W. Schwarz 
[9~T] gave an asymptotic for J2n^x'y( n )~ 1 itself.) 

-[27], 1966- 

De Bruijn establishes Theorem [5] using Rankin's method, a binomial 
lower bound for ty(x,y) and Stirling's formula. 
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-[28], 1963- 

De Bruijn and van Lint show that F(x) ~ G(x) (in the notation of 
[26J), thus answering a question raised by Erdos. 

-[29], 1964- 

De Bruijn and van Lint consider sums of the form Yl n eS(x y) fi n )i w hh 
/ a non- negative multiplicative function (cf. Section [7]). They intro- 
duced a generalization of the Dickman-de Bruijn function (see Section 
El) that around the same time showed up in the analysis of Ankeny and 
Onishi [6] of the Selberg upper bound sieve. 

8.1. De Bruijn's writing style. When de Bruijn introduces a new 
result or method, he also describes to what extent this improves on 
earlier work. Interesting special cases of his result he also likes to 
consider. However, he also discusses limitations and might also describe 
an alternative method. Thus he looks at a topic from various sides and 
gives an honest evaluation of merits and drawbacks of his work. In this 
context, I was struck by the fact that in a lecture he gave at age 90 
[32] he tried to give an honest evaluation of the positive and negative 
sides of his own character! In terms of giving details of proofs he is 
a bit on the brief side (as was usual in mathematical research papers 
written in 1945-1955). Often he promises to come back to some aspect 
in a future paper (and keeps his promise!). 

De Bruijn's book [3D] because of its clarity and pleasant style found 
many adepts. Rather than being theoretical the book proceeds through 
discussing interesting example, several of which are from combinatorics, 
e.g. asymptotics for the Bell numbers^. De Bruijn's book was the 
starting point and inspiration for several later works specializing on 
the use of analytic methods in combinatorics and algorithms, e.g., [TT| 
[811150. 

9. Further reading 

In the first place the reader is advised to consult the book [30] and 
the papers of de Bruijn! As a first introduction to friable numbers I 
can highly recommend Granville's 2008 survey [33]. It has a special 
emphasis on friable numbers and their role in algorithms in compu- 
tational number theory. Mathematically more demanding is the 1993 
survey by Hildebrand and Tenenbaum [64J. The literature up to 1971 
is discussed in extenso in Norton's monograph [83J. For a discussion of 



De Bruijn's book predates the time when the Bell numbers had this name. 
'This paragraph follows closely an e-mail E.R. Canfield sent me. 
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the work prior to 1950 (when de Bruijn entered the scene), see Moree 
[52] . Chapter III. 5 in Tenenbaum's book \§E\ deals with p(u) and ap- 
proximations to *$?(x, y) by the saddle point method, Chapter III. 6 with 
the dual problem of counting integers n ^ x having no prime factors 
^ y. For an introductory account of the saddle point method, see, e.g., 
Tenenbaum [97] . 
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